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Abstract 

Long time ago Pagels and Tomboulis have proposed a model for 
the nonperturbative gluodynamics which in the Abelian sector can be 
reduced to a strongly nonlinear electrodynamics. In the present paper 
we investigate Abelian, static solutions with external charges in that 
model. Nonzero total charge implies that the corresponding field has 
infinite energy due to slow fall off at large distances. For a pair of 
opposite charges the energy is finite - it grows like R a , < a < 1 with 
the distance R between the charges. 
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1 Introduction 



Nontrivial features of QCD such as the asymptotic freedom or the quark 
confinement are due to its gluonic part. Unfortunately, in spite of numerous 
efforts dynamics of non-Abelian gauge fields in a nonperturbative regime is 
not fully understood yet. Popular method of investigating it is the lattice 
approximation. It is the most successful one up to now, yet it has well-known 
drawbacks. There have been several attempts to apply the approach which 
has turned out to be very fruitful in condensed matter physics, that is to 
construct Ginzburg-Landau type effective models which would capture the 
main features of the nonperturbative gluonic dynamics already on the clas- 
sical level. In the condensed matter physics, precisely the Ginzburg-Landau 
models are the main tool for investigating such nonperturbative phenomena 
as, for example, vortices in superfluids and in superconductors, or hedgehogs 
in nematic liquid crystals. 

In the literature there are several proposals for effective models for the 
nonperturbative gluonic dynamics. The most popular ones are the dual su- 
perconductor model [1,2], the colour dielectric model [3,4], and the stochastic 
vacuum model [5,6]. Each of them has its own problems. It seems that the 
satisfactory effective model for the nonperturbative gluonic sector has not 
been found yet, but this does not mean that such model necessarily will be 
very different from the ones already considered in the literature. 

The effective Ginzburg-Landau model is expected to give a good descrip- 
tion of physics already on the classical level (a mean field approximation). 
When constructing such a model, the first step consists in choosing dynamical 
variables, that is the field content of the model. In the case of gluodynamics 
almost all models known to us contain Yang-Mills fields or at least a U(l) 
field if one assumes a kind of Abelian projection, and sometimes other fields. 
Then, it is clear that the effective model for the nonperturbative gluonic dy- 
namics has to be drastically different from typical field theoretical models. 
First, it should not allow for existence of massless, freely propagating vector 
particles (gluons). Second, it should predict absence of coloured states, and a 
confining force between (anti)quarks in colour neutral states. Third, it should 
predict existence of massive glueballs which could propagate freely. Finally, 
the model should explicitly contain a mass scale - for the massive glueballs as 
well as in order to reproduce the trace anomaly of energy-momentum tensor. 
It is a well-known fact that classical Yang-Mills theory with the standard 
Lagrangian obeys none of these requirements. 

In the present paper we investigate a model which has been proposed 
long time ago by Pagels and Tomboulis [7]. In a sense it can be regarded as 
a version of the colour dielectric model without a dilaton type scalar field - 
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the pertinent Lagrangian is a function of Yang-Mills field strength only. The 
model is too simple in order to provide a very good approximation to the 
real-life nonperturbative gluodynamics in QCD, but it looks rather promising 
as far as the four requirements specified above are concerned. Our general 
goal is to study such highly nonstandard field theoretical models in order to 
clarify their physical contents. 

Pagels and Tomboulis calculated electric field of a single point source in 
their model. They found that in comparison with the standard Coulomb 
field it is less singular close to the charge and much stronger far away from 
it. The corresponding energy is infinite due to the large distance behaviour 
rather than the singularity at the origin. This result implies a confinement 
of charges. However, one should also check that a pair of charges Q, —Q 
has finite energy, otherwise also the quark-antiquark pairs would disappear 
from the physical spectrum and the model would be wrong. Because of 
nonlinearity of Gauss law in the model, one can not compute the electric 
field of the pair by a simple superposition. In the present paper, using a 
combination of analytic and numerical methods, we check that the energy of 
the pair Q, — Q is finite if a parameter 5 > 1/4. We also calculate distribution 
of electric field around the charges. It turns out to be very similar to a flux- 
tube, that is it is localised around the straight line connecting the charges, 
but an effective string tension depends on the distance R between the charges. 
The total energy behaves like R a , where a = (45 — 1)/ (45 + 1) and 5 is a 
parameter. Thus, the charges are confined if 5 > 1/4 because then < a < 1. 
The effective string tension behaves like i? a_1 . For 5 = 1/4 we obtain the 
logarithmic behaviour of the energy. 

Interestingly enough, for 5 = 3/4 we obtain the \/R behaviour of the 
energy, which is in agreement with a phenomenological potential found in 
fits to spectra of heavy quarkonia ||. 

Our paper is organised as follows. In Section 2 we recall the model and we 
derive field equations. We also point out that in that model classical vacuum 
is strongly degenerate. In Section 3 we calculate the energy of the system of 
two static, opposite charges. Several remarks about the model are collected 
in Section 4. 

2 The Pagels- Tomboulis model 

Let us denote by JF 2 the standard Yang-Mills invariant^ 

jF 2 = -F^F atlv = B ai B ai — E ai E ai , (1) 
1 We choose the metric in Minkowski space-time with the signature (+1, — 1, — 1, —1). 
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where the field strength F^ u has the usual form 

F; v = d,A a v - d u Al - gf ahc A\A c v (2) 

and E ai = F^, and B ai = —\eikiFu are the non-Abelian electric field and 
magnetic induction field, correspondingly. We will consider only SU{2) gauge 
group, hence a,b,c = 1, 2, 3. The model we would like to investigate has 
Lagrangian C of the form C = ^{T 2 \ where I is a nontrivial function which 
we shall choose shortly. In a connection with QCD models of this kind 
have been proposed long time ago [7], [10], [11]. The Abelian Born-Infeld 
electrodynamics is even older [12], but its aim is to modify only physics of 
strong fields, while in the QCD case the main problem is with weak fields. 
The same remark applies to the quite popular recently non-Abelian Born- 
Infeld actions [13]. 

Pagels and Tomboulis [7] considered Lagrangians of the form 



1 F 2 
2g 2 (t) 



W40 = (3) 



in Euclidean space-time. In formula (|3|) F 2 = F"F®/2 is non-negative due 
to the Euclidean metric, t = In ^f, and g(t) is the effective coupling constant 
determined from the equation 

f m do 

t= 777^ (4) 



where (3(g) is the Gell-Mann, Low function. The Lagrangian (|3|) has been 
proposed in [7] as an Ansatz consistent with renormalization group invari- 
ance of the effective action / rf 4 xL e jj, and also with the trace anomaly. The 
perturbative 1-loop result for the (3{g) function: 

(3(g) = -b g 3 , b > 0, (5) 

where b is a positive constant, gives 

11, 

-2 + 2M, (6) 



9 2 it) gl 



and in consequence 



1 F, 



Leff = ~ + b F 2 \n^. (7) 

Here 1/g 2 is the integration constant, equal to the value of the coupling 
constant at the subtraction point \i. Another case, also considered in [7], 
corresponds to 

P(n) = -Sg, (8) 
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where 5 > is a constant. Now Eq. (M gives 




and 



2d 



J eff 



Notice that 





25 In 


^ 2 y 













l + 2<51n^ + . 



(9) 



(10) 



fill 



hence Lagrangian fll(|) with 5 = b g 2 reduces to (0) when 25 In 4 < 1, that 
is when F 2 is close to the subtraction point //. One may regard Lagrangian 
( [T0|) as a resumation of powers of ln(F 2 //z 4 ) in the effective action functional 
T(Afj) for the Yang-Mills field. From this viewpoint, 5 can be associated with 
the anomalous dimension of the F 2 operator. 

The corresponding effective models in Minkowski space-time are obtained 
by replacing the Euclidean F 2 by T 2 . In order to avoid problems with the 
sign of T 2 and complex- valued Lagrangians we write In in ([?]) as \ In (jjk) , 



and (f) 2<5 as (f )° in ©. 

The model corresponding to L e // (0) was considered by Adler and Piran 
[11] in the case where only non-Abelian electric fields F^ were present - it gave 
a confining force between two point-like opposite charges, at least in Abelian 
sector obtained by assuming that A®{x) = 5^A^(x). However, it turns out 
that if one allows for the presence of magnetic fields, the corresponding energy 
density is not bounded from below. Moreover, Lagrangian (|7|) contains the 
standard kinetic term (c^A"— d v A®) 2 hence the gluons could propagate freely. 
For these reasons, we conclude that the logarithmic model obtained from (0) 
is not satisfactory. 



Let us turn to the model with L e ff given by formula (|10|). When passing 
to Minkowski space-time we change notation a little bit. In Minkowski space- 
time we take Lagrangian of the form 



C 




(12) 



where 5 > 0, instead of Lagrangian (10). The initial coupling constant g has 
been included into A which is regarded as an empirically fixed energy scale 
in the model. 
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Components of energy-momentum tensor corresponding to Lagrangian 
I2T) have the form 



5 

A 8 



2\S 
2 \ 



B 2 + (l + A5)E 2 
T 0l = -(l + 25)e lks E ak B as 



(13) 
(14) 



and 

Tik 



E 2 + (1 + A5)B 2 5 ik -(l + 25) E at l E ak + B ai B ak 



5 



2\« 
2 \ 



(15) 

We see that in this model T 00 > 0. 

Lagrangian ( |i~2"l) does not contain the standard kinetic term for the gauge 
fields. Therefore, it is by no means clear what are physical, propagating ex- 
citations in the model. Let us introduce the non-Abelian dielectric induction 



D° 



— = ( 1 + 25 ) Em {§) 5 



dE° 



(16) 



and the non-Abelian magnetic field: 

dC 



dB Q 



1 + 26)B- (§) 5 



(17) 



The modified Yang-Mills equations following from Lagrangian ([12]) can be 
written in the form 



vy? w = o, v a ob D 



a Ttbk 



a zjbr 



0. 



(18) 



where 



is the covariant derivative. Of course in addition to these equations we have 
the standard non-Abelian Bianchi identities for E a% and B a% : 



V%B U = 0, V a ob B M + s ijk V%E bk = 0. 

These identities are equivalent to formula (2). 
It is clear that any fields such that 



(19) 



(20) 
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obey the field equations because then JF 2 = 0. However, they have van- 
ishing energy-momentum tensor and therefore they should be regarded 
as vacuum fields! The model has enormously degenerate classical vacuum 
state, which includes, in particular, plane waves. In the standard Yang-Mills 
theory, small amplitude classical plane waves correspond to the perturbative 
gluons. It is encouraging that they do not belong to the spectrum of physi- 
cal excitations of Pagels-Tomboulis model. We plan to address the issue of 
physical excitations in Pagels-Tomboulis model in a forthcoming paper [14]. 
In the present paper we show that the model can describe the confinement 
of quarks. 

Formula (16) implies that the dielectric induction is smaller than the 
electric field if < (1 + 28)~ 1 ' , that is the model implies antiscreening for 
weak fields. For strong fields we have the usual screening. 

3 Static solutions with external point charges 

a) The case of non-vanishing total charge 

Equations (18), (19) are even more complicated than the standard Yang- 
Mills equations, which are obtained for 5 = 0. The fact that the classical 
Yang-Mills theory is not confining can be seen already from an analysis of 
its Abelian sector. For this reason, we check the Abelian sector of Pagels- 
Tomboulis model (5 > 0). Let us remind that that sector is constituted 
by gauge potentials A® with only one colour component, identical for all 
fj, — 0,1, 2, 3. For example, we may take 

A; = 5 a 3 A,(x). (21) 

In the Abelian sector, the covariant derivatives in Eqs. (18), (19) reduce 
to the ordinary ones and we may omit the superscript a . Nevertheless, the 
equations remain nonlinear if 5 > - the cases 5 = and 5 > are drastically 
different. In particular, if 5 > 1/4 the Abelian sector of Pagels-Tomboulis 
model is compatible with the confinement of charges. 

Let us consider first a smooth, spherically symmetric distribution jo(r) 
of an external, static charge in a finite region around the origin. The Gauss 
law has the form 

VD= Jo (r), (22) 

and the electric field obeys the condition 

Vx£ = 0, (23) 

with D and E related by formula (16) (we omit the superscript a = 3). 
Because of the nonlinearity, solutions of Eqs. (22), (23) can not be obtained 
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by superposition of solutions for point charges. Nevertheless, spherically 
symmetric solutions follow easily from Gauss law (22). Far away from the 
charge 



Q = / d 3 r Jo (r). 



where 



Thus, both fields vanish when r — > oo, namely \D\ ~ r~ 2 , ~ r - 2 /( 1 + 4<5 ). 
Simple calculation shows that the total energy of the gauge field, / d 3 rT 00 , 
is infinite if 8 > 1/4, due to the behaviour of E at large r. This fact is not 
changed by inclusion of the interaction energy of the external charge with 
the gauge field, / d 3 rj A , because this integral is finite. For comparison, in 
the Yang-Mills case the energy of the gauge field (24) (now with 5 = 0) is 
finite. 

This result suggests that in Pagels-Tomboulis model states with nonvan- 
ishing total charge Q are not physically feasible. In order to turn this sug- 
gestion into a theorem one would have to find solutions of the non-Abelian 
gauss law 

W> K = a»j (r), 

together with the remaining equations (18), (19), without the simplifying 
assumption (21). Even in the Yang-Mills case this is a highly nontrivial task, 
see e. g. jl5] for a review. In particular, Kiskis |16| showed that in the Yang- 
Mills case exist solutions with arbitrarily small (but positive) energy. Such 
solutions in general are not static and have nonvanishing magnetic fields. 
However, his reasoning is based on linearity of Yang-Mills equations in the 
Abelian sector, and therefore it can not be repeated in Pagels-Tomboulis 
model. This is encouraging, nevertheless the question whether finite energy 
solutions with Q ^ are completely excluded remains open. 

b) Dipol-like external charge density 

Effective model for QCD should allow for finite energy colour singlet quark- 
antiquark states 0. Therefore, we shall now consider a dipole-like external 
charge. 

The results of subsection 3a) are valid also when the external charge is 
point-like, jo(r) = Q8(r). Then, the resulting energy density has a singularity 

2 In the case of SU(3) gauge group colour singlets built of three quarks are related to 
the quark-antiquark states because a diquark can be regarded as an antiquark as far as 
colour charge is concerned. 
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at r = which is integrable if 5 > 1/4. Therefore, it is not important for 
our purpose whether the two charges forming the dipole are point-like or 
spatially extended. For simplicity, we assume that they are point-like. Thus, 
we now take 

Jo = qS(x)5{y) [6{z + R/2) - 5{z - R/2)] (25) 

where q > 0, and again we look for solutions of Eqs. (22), (23) in the Abelian 
sector. 

We would like to check that the corresponding total energy £ of the gauge 
field, given by the integral £ = J d 3 rT 00l is finite. First, simple dimensional 
analysis shows that 

, . 86 _ 43-1 

£ = c o| q\ 1+45 A 1 + 4a R 4 *+! , (26) 
where Co is a numerical constant. This formula follows from the fact that 

£ ~ \q\ 1+4,5 A 1 + 4 < 5 , 

as implied by Eq. (22) and formulas (13), (16). The exponent of R is dictated 
by the requirement that £ has the dimension cm' 1 . Thus, it remains to show 
that the constant c is finite. Because singularities at the point-like charges 
are of integrable type, finiteness of Co depends solely on behaviour of the 
fields at the spatial infinity. 



Following Adler and Piran, [T^|, ||T8[, we express the dielectric in- 
duction D by the dual potential C: 

D = V xC. (27) 

Then, at all points where C is sufficiently smooth, automatically V-D = 0. 
However, the presence of source (25) on the r.h.s. of Eq. (22) implies that 
C can not be regular everywhere - it has to have a singularity akin to the 
Dirac string for magnetic monopoles. In our case it is natural to assume that 
the string connects the two point charges. To exploit the axial symmetry 
of the problem it is natural to introduce the cylindrical coordinates (p, (f>, z) 
instead of the Cartesian ones. The coordinates of the point sources are then 



p = 0, z — ±R/2. As in 11171 , we assume that the dual potential C 



can be expressed by a scalar flux function $(p, z), which is defined as follows 



c = ir^ < 28 » 

where is the unit vector tangent to the (j) coordinate line. 
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Equation (23) may be rewritten as 




where 




0, (29) 



(30) 



A 4 




is the dielectric function. The Ansatz (28) reduces Eq.(29) to 

V(<tV$)=0, (31) 

where 

25 

1 \ 1+45 

(32) 

It remains to fix boundary conditions for $. Formulas (27), (28) imply 
that 

$ = for p = 0, \z\ > R/2, 

$ = q for p = 0, \z\ < R/2. (33) 

To obtain, for example, the second line in (33), consider a sphere surrounding 
one of the charges with a small hole around the point at which the sphere is 
punched by the segment of the z-axis connecting the charges. The flux of the 

— * — * 

D field through such surface is equal to line integral of C over the boundary 
of the hole. Shrinking the hole to the point of the intersection we find that 
the value of <3> at that point is equal to the total flux Q of D through the 
sphere. 

We also assume that 

$ -> for p 2 + z 2 -> oo. (34) 

This condition is justified by the expectation that angular dependence of $ 
at the spatial infinity would increase the total energy. Because it has been 
assumed that <!> vanishes on the z-axis if \z\ > R/2, it has to vanish in all 
other directions. 

To summarize, the problem of solving the set of field equations for the 
electric field E with given charge density (25) in three spatial dimensions is 
reduced to one sourceless equation (31) in the region p > 0, z € (oo,— oo), 
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which can be regarded as a cylindrical shell with the outer radius going to 
oo and the inner one to 0. The boundary conditions are given by (33), (34). 

Equation (31) can be rederived from the condition that the flux function 
minimizes the total field energy £ = f d 3 rT 00 under the boundary conditions 
(33), (34). Formulas (13), (27) and (28) give 

1 —2—48 85 2+4(5 

T 00 = -(1 + 45) [2vr(l + 4<5)p]^^ At+^|V$|^, (35) 



where |V$| = y(V$) 2 . One can easily construct examples of the flux func- 
tion which obey the boundary conditions and have finite total energy. For 
example, we may take 

• =%(-, Z + R ' 2 - , : ~ R/2 V (36) 

This function has been obtained from the sum of the dual potentials for two 
Dirac monopoles of opposite charges, and it implies that the Dirac string 
just connects the charges q, — q along the z-axis. The function (36) obeys 
the boundary conditions, but it does not obey Eq. (31) unless 5 = 0. Nev- 
ertheless, the finite value of 8 corresponding to it may be taken as an upper 
bound for the total field energy of the charges q, —q. 

Solutions of Eq.(31) can be found with the help of numerical computa- 
tions. We use the standard iterative procedure [|17]]. Due to the symmetry 
of the problem it is sufficient to restrict our approach to the region z > 0, 
p > 0. The continuous variables p, z are replaced by the computational lat- 
tice with {n p + 1) x {n z + 1) sites; the flux function $ reduces to the discrete 
set of values on the lattice $jj where % = . . . n p and j = . . . n z . The point 
charge q is put on the site of the computational lattice with i = 0, j = n q 
where < n q < n z . The boundary conditions (33), (34) are replaced by: 

$o,j = q, 0<j<n qi (37) 

®o,n q = q/2, 

%,j = 0, n q <j <n z , 

$n p ,j = 0, < j < n z , 

<$>i,n z =0, < i < n p . 

In the first step, sites of the computational lattice are populated by arbitrary 
values. Next, new $jj are computed using a discretized version of the field 
equation (31). The iterative procedure stops when the difference between 
actual value of $jj and the one computed in the previous step is less then 
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a chosen accuracy. By assumption, such approximate solution obeys the 
homogeneous boundary condition (34) already at finite p and z. 

In our computations we used 900 x 900 mesh. Some sample results were 
also obtained for bigger meshes but the results changed insignificantly. We 
have repeated our procedure for some values of the charge separation distance 
R and fixed size of the lattice. The final results were obtained for 5 = 0.75, 
and q = 1.1. The results are presented in Figs. 1-5. Fig. 1 presents the 
flux function $, Fig. 2 the profile of $ for z = 0. From the definition of $ 
the dielectric induction D as well as the energy density e may be derived. 
The rescaled (dimensionless) energy density e = T 00 /A 4 is depicted in Fig. 3. 
Fig. 4 presents the energy density for one of the Coulomb peaks from Fig. 3. 
In Fig. 5 the energy density profiles (e(p) for = z = 0) for various values 
of R are compared. In all figures we use the rescaled spatial coordinates: 
zA — > z, pA — > p and R — > AR.. It is clear that the flux-tube connecting 
the two point charges becomes thicker when the charge separation distance 
increases. We have obtained the same pictures starting from several different 
initial configurations for the iterative procedure. Our numerical procedure 
was also tested for 5 = 0. In this case the analitical solution for <3> is given 
by formula (36). The numerical procedure recovers it. 

The energy density distribution was used to derive the total energy of the 
considered field configuration. As the energy density decreases slowly when 
p 2 + z 2 — > oo it is necessary to collect the energy density from relatively large 
region around the charges. For this reason we have used the lattice with link 
length varying from short near the charges to long for large p 2 + z 2 . The 
energy of configurations for various R was used to check the relation (|26|) - 
if it holds e should depend linearly on y/R for 5 = 0.75. Using numerical 
results for this linear relation one can estimate that c ~ 4.1. 

4 Summary and remarks 

We have investigated the Abelian sector of Pagels-Tomboulis model. If the 
total charge Q is non-vanishing the resulting field has infinite energy, while 
in the case of point-like charges q, —q separated by the distance R the energy 
is finite, proportional to R^A.8 — l)/(45 + 1), provided that 5 > 1/4. 

There are several directions in which one could continue the present work. 
First, we have not investigated stability of our Abelian solutions against 
fluctuations of gauge fields. Because the derivation of screening phenomenon 
given in the standard Yang-Mills case by Kiskis, |16|], does not work if 6 > 0, 
one may hope that a presence of the other components of the gauge field 
{A 1 , A 2 ^) will increase the energy. Also, one may ask about stability against 
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a collapse of the field configuration in the Abelian sector. To answer this 
question one would have to investigate time-dependent solutions. It might 
turn out that additional terms would have to be included in the Lagrangian 
in order to prevent such a collapse, in an analogy to the well-known Skyrme 
term in a mesonic effective Lagrangian. 

Finally, one could use more refined forms of the /3(g) function in Eq. (4). 
The model considered in the present paper should be regarded as the simplest 
one compatible with the asymptotic freedom. 

All problems mentioned above are important, but in our opinion the most 
important and interesting one is to describe physical excitations in the models 
of the Pagels-Tomboulis type. 
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Figure 2: The flux function $ for 5 = 0.75, R = 20 and z = 0. 
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Figure 3: The energy density e for 5 = 0.75 and R = 20. 
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Figure 4: The energy density of the point charge e for 5 = 0.75 and R = 20. 
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Figure 5: The energy density e(p) for S = 0.75 , z = = and R = 10, 20, 30. 
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